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Abstract: 

Comparing the numerically evaluated solution to the leading order GLAP equations with 
its analytical small-x approximation we have found that in the domain covered by a 
large fraction of the HERA data the analytic approximation has to be augmented by the 
formally non-leading term which has been usually neglected. The corrected formula fits 
the data much better and provides a natural explanation of some of the deviations from 
the a scaling observed in the HERA kinematical range. 
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The high-precision data on the deep inelastic structure function F$(x, Q 2 ) coming 
from ZEUS [0J and HI [§] experiments installed in the HERA accelerator have stimulated 
new interest in the properties of the QCD dynamics in this kinematical domain ||. The 
conceptually simplest analysis typically relies on solutions of the GLAP equations which 
are either evaluated numerically or approximated in the small-x region by a suitable 
analytical formula. In particular, Forte and Ball j|, || pointed out that the HERA small- 
x data may be interpreted in terms of the so called double asymptotic scaling phenomenon 
related to the asymptotic behaviour of the GLAP evolution first discovered by DeRujula 
et al.|| many years ago. On the other hand, various groups have been able to fit the 
available data using input characterized by the hard small-x behaviour x~( 1+A ), A > H. 

Notwithstanding the dispute about applicability of the double asymptotic scaling to 
the HERA data, the problem of an analytical approximation to the GLAP evolution is 
interesting in itself ||. A common sense requires that any analytical formula has to follow 
the full solution with an accuracy at least as good as the magnitude of the experimental 
error bars, a goal which may not be so easy to achieve given the high precision of the 
HERA data. Analytical approximations to the GLAP evolution are usually designed to 
be valid in an extreme kinematics of very small x and/or large Q 2 . Starting from this 
observation, we have considered in the present note the question to which extent the 
approximation which leads to the double scaling behaviour can be trusted in the domain 
of moderate Q 2 of the order of ~ 10 GeV 2 which contains a non- negligible fraction of 
the small-x HERA data, even assuming the most favorable case of absolutely soft initial 
conditions. As we shall show in the following, the self-consistent approximation to the 
GLAP evolution in this region indeed seems to require taking into account a formally 
non-leading term which is normally neglected. Alternatively, one may of course try to 
avoid the problem by going to higher and higher Q 2 , which inevitably results in a lower 
statistics, especially in the small x region. On the other hand as we show in the following, 
the improved formula not only provides a much better fit to the small x data, but also 
the appropriately corrected data points are much more consistent with the asymptotic a 
scaling. 

Our analysis is performed at leading order, but its generalization to NLO is not par- 
ticularly difficult. However, before such an extended analysis is performed one has to 
be sure that there are no other sources of corrections which can be even larger than the 
difference between the LO and NLO formula ||. 

At leading order the solution of GLAP equations for moments of the quark singlet 
q s (x,t) parton distribution 
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where t = ln(^-), has the well known form JlO 
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where g(n,t) is the corresponding moment of the gluon distribution and ( = ln(^4^§y). 
The coefficients hx(n) and ^(w) as well as the eigenvalues of the anomalous dimension 
matrix A± arise from the diagonalisation of the flavour-singlet evolution equation. To ob- 
tain the Bjorken-x distribution from (0) one has to take the inverse Mellin transformation 
i.e., evaluate the integral 
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where the integration contour runs to the right to all singularities of q s (n,t). At small-x 
one expects that the bulk of the integral comes from the part of the contour in the vicinity 
of the point n — 0, which can be quantified e.g., by the saddle point approximation. As 
the flavour non-singlet part is subleading at sufficiently small-x values, in this region 
the structure function F 2 (x) should be related to a good accuracy to the singlet quark 
distribution (|) via F 2 (x) = ^q s {n,t) (n f = 4). In Ref.[g § the HERA data were 
analyzed with the help of a further analytical approximation to the integral (f|), based on 
the assumption that in the small- x region and for sufficiently large Q 2 the second term in 
(H), driven by the eigenvalue A_, can be neglected with respect to the first one. In this 
limit the solution of the GLAP equation in the moment space reads 



q s (n,t) = q+(n,t) = ((1 - h 2 {n))q s (n,t ) - hi(n)g(n,t ))exp 



(4) 



The double-scaling solution discussed in Ref. j|, [J is obtained by assuming soft initial 
conditions q s (n,t ) ~ g(n,t ) ~ 1/n, expanding A + (n) and hi(n) (i = 1,2) around the 
point n — 0, retaining only the leading terms, and evaluating the resulting Mellin integral 
(0). As such simple initial conditions typically overestimate the actual magnitude of 
xq(x, Qq) one can introduce an additional parameter xq which has an interpretation related 
to the constant term in the Laurent expansion of the initial conditions around their 
rightmost singularity: 



q(n, t ) = A(Ql) ( - + hi x ) + . . . = A(Q 2 ) - exp (n In x ) + 

\n J n 
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Assuming that xq(x, Q 
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— > 1, with (3 ~ 3 one expects xq ~ e 2 ~ 0.1. 



The reasoning which leads to the approximation (ffl) is based on the series of assump- 
tions which have to be fulfilled by the data under scrutiny. Hence, to be able to justify the 
small-x and large Q 2 approximations we have selected for the further analysis 83 points 
in the domain x < 0.01, Q 2 > 5 GeV 2 from the recently published HI data ||. Although 
for some data in this domain the value of Q 2 is large, relatively many points have much 
lower Q 2 , of the order of 10 GeV 2 . On the other hand, if the perturbative QCD evolution 
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is to be valid, the starting point of the Q 2 evolution should not be much smaller than, 
say, 1 GeV 2 . Then, as the first consistency check we have plotted the distribution of the 
saddle point positions n = (^|) 1//2 , £ = \n(xo/x), for xq = 1 and xq = 0.1, see full and 
shadowed histograms on Figure 1 respectively. It is interesting to see that especially in 
the second case saddle points are not very close to zero, so that one can doubt whether 
the usual small n approximation to the integral (0) is applicable at all. To quantify this 
problem, in the next step we have compared the results of the full solution [|TT[] to the 
evolution equation (|2|) with the double-scaling approximation (|) at Q 2 = 10 GeV 2 , as- 
suming the soft initial conditions for quark and gluon distribution functions q s (n, to) = — , 
g{n, to) — — at Qq = 1 GeV 2 , taking for simplicity A = B. The results are shown on 
Figure 2. The data points correspond to Q 2 = 8.5 GeV 2 , and the normalization of ini- 
tial conditions has been simply adjusted to the data by hand. It turns out that even 
for such simple initial conditions the double- scaling approximation, represented by the 
dashed line, approaches the full solution, depicted by the solid line, at values of x much 
smaller than those experimentally accessible at present. Moreover, in the region covered 
by the data the difference between the full solution and the double- scaling approximation 
is significantly larger than the error bars of the data points, making the applicability of 
the latter to the analysis of the data at least questionable. The situation does not change 



if we assume more realistic initial conditions, say, xq s (x) ~ xg(x, t ) 



-x) 



see Figure 



3 where the appropriate x$ parameter has been taken into account in the double-scaling 
approximation. Our simplifying choice A = B cannot be responsible for the discrepancy 
- a consistent analytical approximation has to be able to accommodate various initial 
conditions without compromising the accuracy. 

Comparing the situation for various values of Q 2 we have found that the situation 
improves considerably when Q 2 increases, or in other words in the above discussion we 
have considered, for the sake of clarity of the argument, a really extreme example. By the 
same token one can suspect that for a reasonable description of relatively low Q 2 data the 
term driven by the eigenvalue A_ should not be neglected. We have tested this hypothesis 
by comparing again the full solution with the approximation in which both the first and 
the second terms in (0) are expanded in n to the leading order, i.e. 



q s (n,t) = ( 0.198 ng s (n, t ) + OAUn g(n,t )) exp 



(6/n- 5.648) C 



+ g s (n,t )exp 
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All real coefficients above result from expansion of corresponding coefficients in @ in n 
and keeping only the first terms. After transformation into x-space the final expression 
for the structure function F 2 {x) to be fitted to the Hl-data in the following reads: 



F 2 {x,t) = — (0.198 A + 0.444 B) exp 
18 
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where 7 = y 12/ Pq. As it can be seen from Figure 4, the formula (|7|) considerably improves 
the quality of the approximation, which in our opinion makes its applicability much better 
justified. 

Motivated by these considerations we have compared fits to the small-x data set based 
on formulae (f|) and fl7|) respectively. For that purpose we have used the minimization 
package MINUIT from the CERN library fll^f . The results are summarized in Tables 1 
and 2. We have kept rif = 4, used Aq C d = 250 MeV and a fixed input scale Ql = lGeV 2 
such that pQCD should be still meaningful. The improved approximation (|7|) results in 
X 2 which is smaller approximately by a factor 3. We note that if we had not restricted Ql 
to a fixed value of Qq = 1 GeV 2 , the fit would drive Ql to much lower values of about 0.5 
GeV 2 , simultaneously pushing xq to values larger than 1, which is certainly disturbing as 
far as the interpretation of Xq is concerned. On the other hand x 2 hardly changes by such 
a shift, and hence we are sure that the values presented in Table 1 represent a sensible 
set of fit parameters. In addition, if the data points in the Q 2 bins around 5 GeV 2 and 
6.5 GeV 2 are removed from the data set, the fit improves once more considerably, see 
corresponding entry in Table 2. 



data range 


0.198 A + 0.444 B 


Xq 


X 2 /d.o.f. 


Q 2 > 5 GeV 2 , x < 0.01 


1.60 ± 0.06 


0.14 ± 0.01 


1.57 


Q 2 > 8.5 GeV 2 , x < 0.01 


1.41 ± 0.06 


0.19 ± 0.02 


1.35 



Table 1: LO fit with double-scaling solution 



data range 


A 


B 


Xq 


X 2 /d.o.f. 


Q 2 > 5 GeV 2 , x < 0.01 


0.83 ± 0.09 


2.47 ± 0.18 


0.15 ± 0.02 


0.49 


Q 2 > 8.5 GeV 2 , x < 0.01 


0.80 ± 0.10 


2.30 ± 0.18 


0.18 ± 0.03 


0.35 



Table 2: LO fit with improved solution 



What is the influence of the correction taken into account in (|7|) on the double asymp- 
totic scaling considered in @, [J? To see this we have compared the LO double asymptotic 
scaling formula which results from @ with the data in the domain x < Xq, Q 2 > 8.5 GeV 2 
before and after the correction due to the extra term in (^) is taken into account. The so 
called cr-scaling predicts Q that after rescaling the F 2 -data by the factor 

R F = iV F p v /^exp(-2 7o- + 5 + a/p) (8) 

where 6 + = (11 + J^n/)//?o, one should observe constant behaviour at large p, where 
P — \/s/C ) a — vTC; 7 — \/12/A>) an d Np is an arbitrary normalization constant. The 
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results are presented on Figure 5. The upper plot corresponds to the original double- 
scaling formula (§). For the lower plot our fit result for the non-scaling contribution 
arising from the formally non-leading term in (0) has been subtracted from each data 
point before the rescaling. For p > 1.5, points on the lower plot clearly have much smaller 
spread around the horizontal line. The p scaling predicts that on the plot 1ii(Rf> F%) 
versus a the data points, with the rescaling factor given by 



Rpi = N F ip^fa exp(<5 + aj p) 



(9) 



should lie on the straight line with the slope equal to 27. To see the effect of ([/]) on the 
p scaling hypothesis we have displayed on Figure 6 the scaling behaviour before (upper 
plot) and after (lower plot) the non-leading correction has been subtracted from the data 
points. Although both plots clearly differ in the domain of a < 1, the scaling behaviour 
is hardly influenced for o > 1.5. The dashed line represents the expected behaviour 
according to the first, scaling term in (fil) i.e., 



(10) 



R F ,F 2 = Ny/chp'ya) 

with the constant iV adjusted to the normalization of the large a data points. 

For large a the data plotted in Fig. 6 should follow a straight line with slope 27 [JL3 
Restricting our data range to Q 2 > 8.5 GeV 2 , x < 0.01 we have determined, using 
MINUIT as well, the slope before and after the non-leading correction has been subtracted 
from the data points. Our results, summarized in Table 3, clearly show that with the non- 
leading correction taken into account we get better agreement with the value expected 
from the theory, 2jth = 2.4 





data range 


2 7 


double scaling solution 


a > 1.0 


2.15 ± 0.04 


improved solution 


a > 1.0 


2.47 ± 0.05 



Table 3: Slope 27 obtained from a fit to i^-data 



To summarize, starting from the observation that the experimental errors set the 
magnitude for the minimal sensible accuracy of an analytical approximation to the GLAP 
evolution, we have found that the description in the region of Q 2 ~ 10 GeV 2 improves 
considerably after a formally non-leading, and thus usually neglected term is included. 
The improved formula provides a better description of the available small x HERA data 
and it improves the observed a scaling. 
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Fig. 1 Distribution of saddle points n = (^|)^ 2 over the analyzed HERA data points 
[[|. The shadowed and filled histograms correspond to Xq = 0.1 and Xq = 1, respec- 
tively. 
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2 Comparison of the full evolution (solid line), equation (§) with the approximate 
double-scaling solution (|]) for the data around Q 2 = 10 GeV 2 [0] , and the soft initial 
conditions q(n,t ) = g(n,t ) = 0.45/n at Ql = 1 GeV 2 . 
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Fig. 5 The a scaling formula from Ref.[Q, |5| versus the small x HERA data 0. For the 
lower plot the formally non-leading contribution taken into account in (0) has been 
subtracted from each data point before rescaling of the data M according to (pf). 
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Fig. 6 The p scaling formula from Ref . 0, versus the small x HERA data [Q] . For the 
lower plot the formally non-leading contribution taken into account in (|7]) has been 
subtracted from each data point before rescaling of the data @ according to @. 
The dashed line represents the expected behaviour according to the first, scaling 
term in (0). 
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